Orthogonality catastrophe and decoherence in a trapped-Fermion environment 
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The Fermi edge singularity and the Anderson orthogonality catastrophe describe the universal physics which 
occurs when a localized scattering potential is suddenly introduced in a Fermi sea leading to a brutal disturbance 
of the ground state. We demonstrate that the effect can be seen in the controllable domain of ultracold trapped 
gases by providing an analytic description the out-of-equilibrium response to an atomic impurity, both at zero 
and at finite temperature. Furthermore, we link the transient behavior of the gas to the decoherence of the 
impurity, and, in particular to the amount of non-markovianity present in its dynamics. 
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A Fermi gas may be shaken-up by the switching of even 
a single, weakly interacting impurity, producing a complete 
rearrangement of the many body wave-function that, as a con- 
sequence, loses essentially any overlap with the initial, unper- 
turbed one. This is the essence of Anderson's orthogonality 
catastrophe [ 1|, witnessed by the singular (edge-like) behavior 
of the energy distribution of the impurity induced excitations. 
An example of how such a many-body effect comes into play 
is provided by X-ray photoemission spectra from most sim- 
ple metals [2, 3 1, where the expected sharp symmetric peak at 
the binding energy of a core level is converted into a power 
law singularity. Similar patterns have been observed in elec- 
tron emission via X-ray absorption and Auger neutralization 
from carbon based nanomaterials (|4l, and recently predicted 
for quantum dots [ 5 ] . Fermi edge resonance and orthogonality 
catastrophe have been also revealed by non-equilibrium cur- 
rent fluctuations (shot noise) in nano-scale conductors O, and 
enter prominently the physics of phenomena as diverse as the 
Kondo effect [7 1 and the scattering or sticking of a low-energy 
atom or ion on a metal surface ||8] [9] . 

Recently, it has been proposed to observe this universal 
physics in controllable ultracold atomic setups where the sin- 
gular behavior may be probed either in the time domain 
by Ramsey interference type experiments, performed on the 
impurity atom [10], or in the frequency domain by radio- 
frequency spectroscopy ifTTl . However, an analytic frame- 
work for the case of a trapped Fermi gas is lacking. In 
this letter we provide such an analytic description, and dis- 
cuss the transient response of a harmonically trapped Fermi 
gas following the 'sudden' switching of an embedded two- 
level atom excited by a fast pulse. The switching of the lo- 
cal impurity produces a local quench of the gas, giving rise 
to the Anderson catastrophe. We propose a study of the 
Fermi-edge physics at zero and finite temperature and both 
in the frequency domain, by looking at the excitation spec- 
trum of the gas, and in the time domain, by studying the 
decoherence dynamics of the impurity. This allows us to 
link the Fermi edge behavior of the excitation energy dis- 
tribution to the dynamics of the impurity, and, in particular, 
to the Loschmidt echo 1T21 [131 . Furthermore, we investi- 



gate the non-Markovianity of the induced dynamics using re- 
cently developed tools I14H17I1 . which have been employed to 
study open systems in different environments, ranging from 
spins lfl8l to Bose-Einstein condensates [ 19 1, and experimen- 
tally tested in optical set-ups [ 20 21|. We find that the non- 
Markovianity of the decoherent dynamics provides a novel 
interpretation of the essential physics of the process. Specifi- 
cally, we consider a gas of non-interacting cold fermions con- 
fined by a one-dimensional trapping harmonic potential of fre- 
quency lj, described by the Hamiltonian Hq = J2 n EnC^Cn, 
with C n being the annihilation operator for the n-th single 
particle level of energy e n = hui(n + 1/2). An impurity is 
present in the gas, constituted by a two-level atom with en- 
ergy separation hfl, localized at the minimum of the potential 
well. The fermions lie in their equilibrium configuration, set 
by Ho, until the impurity is excited and a sudden perturba- 
tion V(t) = V9(t) is felt by the gas. To mimic a very strong 
difference in scattering length depending on the internal state 
of the impurity, we model it by a spatially localized (delta- 
like) potential, activated by the population of the excited state, 
V = V(x) (g) |e) (e|, where V(x) — ttVoXoS(x). Here, Vq is 
the coupling constant and xo is the characteristic length of the 
harmonic potential. Only even one-fermion states (with n = 
2r) are involved in the shake-up process, being coupled by the 
matrix elements V„> = V (-l) r+r ' 2 T 1 (r) 1/2 2 Y 1 (r') 1/2 , 
with denoting the Euler's gamma-function ratio 2^1 (r) — 
r (r + 1/2) /r(r + 1). Without loss of generality, we may as- 
sume the Fermi number to be even so that the Fermi energy 
reads e F = Hlj (2r F + 1/2), with level occupation numbers 
given by the Fermi factors /,. = (l + e 2/3huir ^ 

A key quantity for the following is the so-called vacuum 
persistence amplitude 

vp(t>0) = (e^ ot e-^ o+9 >), (1) 

with (...) denoting the thermal equilibrium average over the 
unperturbed environment states, v$ is the probability ampli- 
tude that the environment will retrieve its equilibrium state at 
time t, after the switching on of the perturbation and, as dis- 
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cussed below, it gives the decoherence factor for the impurity. 
The Fourier transform 



a = 0.4 a = 2 



ME) = f 

J — I 



dt 



h Vf){t), 



(2) 



subject to the analytical constraint vp (t < 0) = v% (—t), is the 
excitation spectrum relative to the unperturbed energy. 

In the interaction picture we get 

vp{t) = (Te-s/o dt '^')\ v{t) = e 4% e 4 A »', 



(3) 

which, by virtue of the linked cluster theorem, reduces to an 
exponential sum of connected Feynmann diagrams, vp (t) = 
exp{A /3 (t)}, with: 
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After introducing the unperturbed propagator 

ihG?{t) = e-^ 1 / 2 ^ [9{t) / rF _ r - 9{-t) / r _ rF ] , (4) 

and applying Wick's theorem, we focus on the loops with one 
and two vertices, namely, (f) and A^i). These diagrams 
contain separable products of arrows (G@) connected to ver- 
tices (V rr i). In terms of the auxiliary functions 



we may express the diagrams as 

\P(t) = ~aV X+(0), 



(5) 



(6) 



2 ft ft 1 



h 2 Jo 



dt 1 I dt"\1(t") \t(-t"). (7) 



We will see in the following that the relevant energy scale 
for the problem is set by ^/hu>s f , involving both the level 
separation and the Fermi energy. Therefore, we introduce 

aV 2 

a = 2 h J' e as a sensible parameter, ranging typically from 
to 1 (a accounts for the spin degeneracy). The contribution (|6]l 
is then written as Af (t) = —itE^/h, where 



E{ = ^2ahuje e aJ2 2^1 (r) f r 



r=0 



is nothing but the first-order correction, as obtained by using 
Rayleigh-Schrodinger perturbation theory, to the unperturbed 
energy Ejj = J2 n ?n f n /2-r r > which is plotted in Fig. 
Then, A^ (t) brings a phase factor to vp (t), which corresponds 
to shifting the spectrum vp(E) by E^, whose behavior vs e F 
is shown in Fig. |T^5. We notice that e\ is 1 to 3 orders of mag- 
nitude larger than £jf , depending on j3 and e F , while temper- 
ature plays an appreciable role in both Eq and E± for e^/huj 
below ~ 100, and /3Hlj less than ~ 0.1. 
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FIG. 1. (Color online) Unperturbed equilibrium energy Ep of a spin- 
1/2 gas into a harmonic trap (panel A), and perturbation corrections 
Eg (panel B), Eg (panel C), due to the sudden switching potential 
V(x), Energies are reported in units of huj vs e ¥ /hu) for values of 
Phw in the range of 0.001 to 1, and in the fihuj — > oo limit (absolute 
zero). Finite temperature summations are computed from Eqs. <[8j 
and 1 12 1, with a = 0.4, using a high energy cut-off r < 10 2 r F . 



Performing the time-ordered integrals in Q, we rewrite the 
two-vertex connected graph as 



in which 



'■pr 



E 



2lu 



Jr—rrfr, 



iTi (r) 2 ri (/) 



(9) 



(10) 



1 — 2i(r-r')tu 
Ar'(t) = 2 Tl (r) 2 T 1 (/) — -g—. (11) 

(r — r ) 

Here, we may separate the sums over even-state labels (r,r') 
into three contributions with well defined trends, h^(t) — 
A| s (i) + A^ G (t) + A^ p (i), where the subscript stand for Shift, 
Gaussian and Periodic, respectively. In particular, the off- 
diagonal summands in (|9]l that multiply (p rr i provide a (fur- 
ther) shift to the environment energy, the second-order cor- 
rection Af s (t) = —itE$ /%, with 



^ ^ ^rr' fr—r f fr F - 



(12) 



r,r'— 
r^r' 



In Fig.[T|c, we notice that E^ takes absolute values similar to 



(8) Ei, but it is more sensitive to temperature, for (3tua < 0.05. 
The diagonal elements of Eq. Q yield a quadratic power law 



A 2gW = -a^t 2 gp, 



(13) 



which produces a gaussian damping in vp(t), and a gaussian 
broadening in i/p(E), depending on the coefficient 

oo 

90 = E ( r ) 2 ^- r F ^-f ^ 
r=0 

The broadening parameter is weakly influenced by the 
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Fermi energy, but strongly affected by temperature, inducing 
changes of various orders of magnitudes for (3fku < 2 and 
tending to atv 2 in the zero temperature limit (Figj2]\). 

The most important content of the second diagram (t), 
giving a non trivial structure to vp(t), arises from the series 



2hw 



OO 

E 

r,r'— 



(15) 



being periodic in time with frequency uj (see Fig. |2jJ). The 
zeroes of this sub-graph, when combined with the gaussian 
damping ( p"3j ), yield modulations in the vacuum persistence 
amplitude which, as discussed below, are a signature of non- 
markovianity. 
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FIG. 2. (color online): (A) Gaussian broadening of eq. 
units vs Pfrw for r F = 5. 50, 100, and a = 0.1. (B) Modulus of the 
sub-diagram A§ E (t) |l5| , displaying periodic behavior with respect 
to tot, represented for f3hu> = 0.01 — 1.5, r F = 100, and a = 0.1. 



The Gaussian power law of Eq. ( 13 i and the periodic con 



tribution of Eq. ( 15 1 can be approximated as 



<(*) 
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1 



^ (-l) m m 

/ j e 2m/3 _ ^ 



(16) 



Af P (*) = In 



e 2r °" - 1 

e 2uj(T +it) _ ^ 



E ln 

m— 

oo 

E 1 * 



l _ e 2(it-(m+l)f3h)u> 
e 2(it+m/3h)uj _ ^ 



(17) 



where To is a regularization parameter (needed since A^i) 
is singular at tq —> + ), that we can interpret as the typical 
time-scale over which environment transitions occur. Notice 
that thermal fluctuations introduce other characteristic times 

T m = m(3h. 

Numerical calculations of vp(t) at Fermi numbers below r F ~ 
100 may be run straightforwardly, using Eqs. ( 13 1 and ( 15 1, 



but an analytic approximation is easily obtained as well: leav- 



ing aside the shifts due to Eqs. d8b and ( 12 1, the non trivial part 



of the vacuum persistence amplitude is 



(18) 



Of particular interest is the behavior of 1 1^/3 (i) | exhibiting 
spikes at cut ~ n, 2ir, . . ., which become more and more pro- 
nounced with increasing fihu), see left panels in Fig. [5] 
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FIG. 3. (color online): Absolute value of the decoherence factor 
(left panels, A,C,E) and absorption spectrum v'p(E) (right 
panels, B,D,F), calculated from eqJTSl by numerically computing 
the sub-diagrams of eqs. < | 1 3| > and (15), for f3tho = 0.1 — 10, r F = 
100, and ol = 0.1 - 1.05. 



The periodicity in the time domain is reflected in the absorp- 
tion spectrum Vp{E) that offers an asymmetric, broadened, 
signature of the singular behavior of the fermion environment. 
The monotonic structure turns into a sequence of sub-peaks, 
separated by 2huj and related to even-level transitions in the 
gas as f3huj gets above ~ 0.5 (see Fig. |3jJ). These features are 
observed for any r F in the range of 5 to 100 |22|. Eqs. ( 16 1 
and 



17 1 give the following form to the vacuum persistence 



amplitude at zero temperature, 



^-at 2 bj 2 /2 



e 2ui(T +it) _ ^ 



(19) 



which we checked to accurately reproduce the exact numerical 
result for Fermi numbers larger than 10, for suitable values of 
the cut-off parameter, say, t lu ~ 0.01 (see Fig. |4|V). 

We now compare our results to the one-dimensional free- 
fermion case. To this end, we keep a and e F constant, letting 
hid — >• with 2r F « e F /hu> being large but not infinite. No 
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gaussian damping occurs in this case, while (t) tends to 

1 



AmndO) = log 



U/t q + 1 



(20) 



yielding the Mahan-Nozieres-De Dominicis (MND) propaga- 
tor ^mnd(0 = (it/ tq + l)~ a , originally calculated for a sud- 
denly switched core-hole in a free electron gas [2, 3J. The 
parameter fi/r coincides with the effective energy range of 
the core-hole potential around the Fermi energy, and a is 
the MND critical parameter. Interestingly, in our case the 
a-parameter bears the same range of criticality as the MND 
parameter (a — — 1). Finally we notice that temperature 
effects are mainly contained in 



n 



I _ e -2(m+l)/3hu> 
I _ e 2(it-(m+l)f3h)uj 



; 2w(T +it) _ X 



1 



D 2(it+m/3h)u) _ X 



which provides an excellent approximation of Eq. ( 19 1 



Up to now, we treated the response of the Fermi gas with- 
out any reference to the dynamics of the impurity, that has 
just been assumed in the excited states for t > 0. If, instead, 
the two-level atom is subject to (say) a ir/2 pulse and pre- 
pared in the superposition (\g) + |e))/v2, it experiences a 
purely dephasing dynamics such that its state at later times is 
Pmp(t) = (Iff) (<?| + |e) (e| +»//,(*) Iff) (e|+h.c.)/2. Thedeco- 
herence factor that enters the off-diagonal elements is just the 
persistence amplitude that we obtained before, going to zero 
at long times due to the Orthogonality Catastrophe. In the 
theory of open systems, one typically uses a related function, 
the so called Loschmidt echo L(t) = \vp(t)\ , which gives 
a measure of the environmental response to the perturbation 
induced by the system |[T2l [131 123| and which, as shown in 
Ref. fPTl . is linked with non-Markovianity of the open system 
dynamics. The amount of non-Markovianity of a dynamical 
map can be evaluated in different manners lfl4"HT6l . which are, 
however, essentially equivalent for a purely dephasing quan- 
tum channel lfT7ll24l . By adopting the definition in terms of 
information flow given in [14|, one finds 



(21) 



Af = > \VL(t U AX,n)\ ~ |VL(W,n)| 



where the sum is performed over all maxima and minima of 
the Loschmidt echo. Using our previous results for the am- 
plitude vp(t), we can then obtain the non-Markovianity of the 
dynamics of a two-level system in a trapped Fermi gas. The 
results are shown in Fig. [4ji, where we see that Af depends 
on the temperature and on the critical parameter a. In par- 



ticular, it goes to zero both for large temperatures (due to the 
fact that thermal fluctuations suppress oscillations in the per- 
sistence amplitude) and for a > 1. 

In the latter case, excitations are generated at every energy 
scale in the fermion gas, as witnessed by the fact that the 
spectrum becomes structure-less. This implies that the gas 
becomes more and more stiff (in the sense that it is not able 
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FIG. 4. (color online): (A) Absorption spectrum i/'p(E), calculated 
numerically from eq. l |18[ with fffiw = 0.1 — oo, r F = 5, and a = 
0.1, and analitycal approximation u'^^E) obtained from eq. | |19| ; 
(B) Non-Markovianity measure as a function of the critical parameter 
a for various temperatures. 



to react on the impurity any more). This explains why Af is 
zero: the open system does not receive information back, its 
Loschmidt echo decays monotonously and thus the dynamics 
is Markovian. As a result, we conclude that a non-Markovian 
dynamics can be characterized, in our case, by the appearance 
of specific spectral features in the excitation energy distribu- 
tion |25l . 

We conclude with two final remarks. The first one is that 
the spectral distribution of energy excitations we have ob- 
tained above coincides with the so called work distribution 
function, which is a central quantity in non-equilibrium pro- 
cesses [26 1 . In the set-up that we have described above, it is 
simple enough to imagine a 'reverse' protocol, with the fermi 
gas brought to thermal equilibrium in the presence of the im- 
purity (i.e. with the two-level atom in the excited state) which 
is then switched off. The comparison of the work distribution 
functions in the direct and reverse protocols would lead to a 
direct experimental test of the Crooks relation in the quantum 
regime (27) . Finally, we remark on the experimental realiza- 
tion of the model we have described. Many experiments have 
recently dealt with the effects of impurities in trapped Fermi 
gas (28 1, and state-dependent scattering lengths have been dis- 
cussed 1 29 ]. This would lead to a direct test of our theory. 
Another viable candidate could be a gas of hard-core bosons 
in one-dimension, where the Loschmidt echo is equivalent to 
that of the corresponding Fermi gas |30| and in which impu- 
rities have recently been experimentally generated ll3"Tl . 
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Supplementary Material 



APPENDIX A: IMPURITY POTENTIAL 

The impurity, <5-like potential V(x) employed in this work 
is localized at the minimum of the harmonic well, thus, it in- 
duces excitations that connect only unperturbed one-fermion 
states labelled by even numbers n = 2r. The diagonal 
matrix elements of V(x), shown in Fig. |5jV, entirely deter- 
mine the one-vertex graph of Eq. (6) of the main text, rep- 
resenting first order corrections to the single-particle ener- 
gies e n = TvjJ (n + 1/2). In addition, both diagonal and off- 
diagonal matrix elements of V(x), the latter shown in Fig.Bfe, 
appear as absolute squares in the two-vertex graph A^yt). 
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FIG. 5. (color online): (A): First order energy shifts for the occupied 
fermion levels e n of the trap. The continuous limit representation of 
V nn is obtained by letting n take any real value between and 2r F , 
which corresponds to the cu — > 0-limit. (B): Second order couplings 
between even one-fermion states that enter the two-vertex graph of 
Eq (7). 



APPENDIX B: NUMERICAL COMPUTATIONS 

The quantities involved in our model for the vacuum per- 
sistence amplitude Up(t) are: 

• the energy of the fermion gas in the harmonic trap (Eq), 
together with the first and second order shifts due to the 
impurity potential (e.g., Ef and E^ , reported in Eqs. 
(8) and (12) of the main text); 

• the gaussian damping coefficient gp, entering the sub- 
diagram Ajg, given as Eq.(14) and (13) in the main text, 
respectively; 

• the Fermi edge loop A^ P reported in Eq. (15). 

These contributions contain summations running over all one- 
fermion eigenstates of the trap that we truncated by us- 
ing a high energy cut-off e CUT ^> e F . To fix e CUT = 



hui (2r CU T + 1/2), we performed convergency tests by chang- 
ing r CUT in order to have a maximum instability error below 
0.1%. We observed that accurate estimations of E^ and gp re- 
quire values of r cux in the range of 10 3 to 5 x 10 3 for r F below 
~ 100 and /3Hu) larger than ~ 0.1. In particular, the plots of 
Fig. 1C and Fig. 2 A were generated with a energy cut-off of 
25 - 500 e F . As for E%, E{, and A^ p , we found out a cut-off 
of ~ 10 e ¥ to be sufficient in the investigated range of fermion 
numbers and temperatures. Accordingly, in fig. 1A, fig. IB, 
and fig. 2C, we included up to 10 3 one-fermion states. 



APPENDIX C: ANALYTICAL APPROXIMATION FOR THE 
TWO-VERTEX GRAPH 



The expression of the free fermion propagator, reported in 
Eq. (4) of the main text, and the power series expansion 



L = { 1/2 " 



r < 

r = (22) 
r > 



of the statistical Fermi factors f r -r T , allow us to manipulate 
the connected diagrams entering vp(t), and find analytical ex- 
pressions for h-2 Q (t) and K 2v (t). 

In particular, we may rewrite the gaussian damping coefficient 
introduced in Eq. (14) of the main text as 



9H 



2F1 (r F ) , „, ; 



+ > (-l) m mg^(r F ), (23) 



m — 1 



in which the finite temperature summand 



f m (r F ) = 2 r x (r F ) 2 3 F 2 (r F ,e 2m0 " n ) 

- 2T1 (r F + l) 2 e - 20muh 3 F 2 (r e + 1, e ~ 2m ^ h ) 



-2 K (e 2ml3ujh ) 
depends on the generalized hypergeometric function 



(24) 



3 F 2 (r F ,z) = 3 F 2 (l,r F + l/2,r F + l/2|r F + l,r F + l|z), 

(25) 

and the complete elliptic integral of the first kind K. 

We consider a many fermion environment (r F 3> 1), approxi- 
mating gf^ with the large-r F expansion 



1 

4r F 



E 



(-l) m m 

e 2m/3 _ I 



(26) 



that works extremely well for r F > 10 and f3 > 0.5. The 
asymptotic form of gp leads to the result reported as Eq. (16) 
in the main text. 

As for the Fermi-edge component Agp, we work on the aux- 
iliary functions introduced in Eq (5) of the main text and use 



7 



the power series expansion reported in Eq.(22 1 above, to write 



xi(t) 



m=0 

oo 



r F -l 
r=0 



± e ±2writ 



2 ri (r F ) 



m=0 



(-l) r 



/ J V-m-l ' 
r— r F + l 



(27) 



,.,„ - 2 r\ (r) e^Mr-M, and l~ m = l+ m+v Next, 



with Z+ 
we consider the summations 



r F -l 

E » r i w 

r=0 



z rF 2 A(r F ,z), 



and 



E - 

r— r F +l 



r x (r) z r = ^ +1 2 A(r F + M), 



(28) 



(29) 



holding for any z ^ 1, where is the regularized 2-F1- 
Hypergeometric function 

2 F 1 (r F ,z) = 2 r x (r F ) 2 F 1 (l,r F + l/2,r F + l;z), (30) 

having the asymptotic behavior 

„-V2 



2 Fi (r F > l,z) 



1 



°( r " /2 ) 



(31) 



With the representation ( |27| i, we need to add an imaginary 
time regularization to the t" -integral in the two-loop term, 
e.g., we have to shift the t" integration domain by itq to pre- 
vent 2 F 1 (r F ,e ±2urlt ") and 2 Fx(r F + l,e ±2cjrit ") from be- 



ing singular. We, then, insert eq. ( 27 1 in the expression for 
A2 (t), compute the r-summations, and use the large r F expan- 
sion |3~T| i. Finally, we perform the t' and t"-integrals, exclud- 
ing terms proportional to t and i 2 . What is left is a combina- 
tion of logarithmic and polylogarithmic functions, dominated 
by the Fermi-edge term reported in the Eq. (17) of the main 
text. Indeed, as shown in Fig. [6] the zero temperature form 



A£(t)«ln 



1 



(32) 



obtained from this procedure is in excellent agreement with 
the numerical calculations reported in Fig. 2B for r F ~ 100 
and f3hui larger than ~ 1. The correctness of such an approx- 
imation is also attested by the comparison of Fig 5A. 



APPENDIX D: DECOHERENCE FACTOR AND 
ABSORPTION SPECTRUM 

We performed numerical calculations of the decoherence 
factor and the environment absorption spectrum 

vp(E), by selecting different fermion numbers (r F = 5 — 
100), coupling parameters (a = 0—1.1), and thermal en- 
ergies (fihw = 0.001 — 10). In the main text, we have pre- 
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FIG. 6. (color online) Absolute value (A) and phase (B) of the Fermi 
edge diagram Ajp, for a = 0.1 and 2r F = 100. Numerical com- 
putations from Eq. (15), for /3foj = 0.9, 1.5 are compared with the 
ft — > 00 approximation obtained from Eq. l|32[>. 




FIG. 7. (color online): Absolute value of the decoherence factor 
(left panels, A,C,E) and absorption spectrum v'p(E) (right 
panels, B,D,F), calculated from Eq. (18) of the main text by nu- 
merically computing the sub-diagrams of Eqs. (13) and (15), for 
Phcu = 0.1 - 10, r F = 5, and a = 0.1 - 1.05. 



sented an application to a Fermi gas of 200 particles (in Fig. 3) 
where the Fermi-edge behavior, superimposed on a gaussian 
damping trend, appears as a sequence of spikes in and 
an asymmetric peak structure in vp(E). Such features be- 
come more and more marked with decreasing temperature, 
which reduces the effect of the gaussian coefficient gp (see 
Fig. 2A), whereas they disappear for a > 1. Similar con- 
siderations hold for environments containing low numbers of 
fermions (see Fig. [7]) in which shake-up effects are even more 
visible because of the decreasing of gp with decreasing e F . 
On the other hand, the approximation given by Eq. (18) in the 
main text works extremely well with environments containing 
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large fermion numbers. Nevertheless, we noticed that values 
of r F larger than ~ 10 are sufficiently well reproduced by the 
analytical model, at any temperatures. 



A «=0.1 B 




FIG. 8. (color online): Absolute value of the decoherence factor 
1 1^9 (t) | (A) and absorption spectrum v'p(E) (B), calculated from 
Eq. ( 1 8) of the main text by numerically computing the sub-diagrams 
of Eqs. (13) and (15), by including 2000 m-products, for /3hui = 
0.1 - 1, r F = 50, and a = 0.4. 



